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Abstract 



A method is described for estimating muon track momentum from 
the distribution of hits along tracks in dense calorimeters due to multi- 
| pie Coulomb scattering (MCS). The method requires only well-aligned 

■ conventional drift chambers or similar tracking devices and can be 

0^ | implemented with or without a magnetic field present. Neutrino ex- 

^ ■ periments can use this technique to increase acceptance for large angle 

1} \ charged-current scattering events. Resolutions in a typical detector 

geometry vary from ~ 10% at p = 20 GeV/c to ~ 50% for p = 500 
QJ ' GeV/c, if the number of hits on the track is sufficient. 
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1 Introduction 

Neutrino detectors are commonly constructed as long dense calorimeters to 
maximize interaction rate. This geometry results in loss of acceptance for 
charged current interactions from muons that exit the sides before they 
reach the spectrometer that is typically immediately downstream of the 
calorimeter. In toroidally magnetized calorimeters, muons can exit before 
sufficient BdL is accumulated to measure momentum, or they may leave a 
large fraction of their track length in the central hole of the toroid. Losses 
are greatest at high values of Bjorken scaling variable x and inelasticity y. 

Fortunately, these detectors are often instrumented with a large num- 
ber of tracking chambers to determined the neutrino interaction vertex and 
muon scattering angle. Resolution on this angle can be dominated by mul- 
tiple Coulomb scattering (MCS) up to TeV energies. Strong dependence 
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of MCS error contribution on momentum and the large number of hits on 
a track in neutrino detectors permits a different momentum determination 
scheme. The procedure dates from the late thirties Q, has been used in 
many emulsion experiments and is still used in balloon-borne cosmic ray 
experiments with a variety of tracking technologies |||. It entails a straight 
line fit to a muon track that varies slope, intercept, and momentum such 
that the probability distribution for the observed pattern of hits is maxi- 
mized. The MCS-based momentum estimation does not require a magnetic 
field and allows for a substantial recovery of the acceptance loss from exiting 
muons. Reasonable resolution can be obtained for muons with momenta up 
to several hundred GeV/c using a straightforward track finding and fitting 
algorithm 

The following sections describe the procedure in more detail and the re- 
sults of calculations for a detector geometry consisting of N identical track- 
ing chambers with spatial resolution do separated from each other by a 
constant thickness A of material with radiation length Xq. The calculations 
are tested with a GeantQ] Monte Carlo simulation of the NuTeV neutrino 
experiment!! at Fermi National Accelerator Laboratory. This experiment, 
chosen for its "typical" neutrino detector, is briefly described in Appendix 
[A|. It has parameter values of N < 42, do = 0.05 cm, A = 42.4 cm, and 
Xq = 3.45 cm for the purposes of this paper. A forthcoming publication will 
provide results of application of the procedure to NuTeV data. 

2 Tracks in a Dense Detector 
2.1 x 2 Based Momentum Estimation 

Consider fitting a small-angle muon track to a straight line in a dense 
calorimeter instrumented with many equally spaced tracking detectors (as- 
sumed to be drift chambers for the sake of discussion). This may be accom- 
plished by minimizing a x 2 function that compares measured hit positions 
to a linear trajectory, 

X 2 = (y~ %zi - yozo) V _1 (p) (y- Q zi - y z ) , (1) 

with respect to the slope 0$ and intercept t/o- Here, y and z\ are the ./V mea- 
sured y, z points, and zq is an N dimensional vector with all of its elements 
equal to unity. The covariance matrix V(p) contains constant contributions 
from chamber resolution and momentum dependent terms from multiple 
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Coulomb scattering (MCS): 

Vij = o%8ij + Sij(p), (2) 
where the scattering matrix element, usually attributed to Fermi 0], is 



min(ij) 2 

sM= E ft 

k=l P k 



A 2 A 

k + -^-Oi - 2 fc + Zj - Z k ) + (Zi - - z k ) 



3 2 

In these expressions for the covariance matrix, a® is the drift chamber resolu- 
tion, A/% is the distance in z between hits k and k— 1, Zi is the distance from 
the track start to the i th hit, and p k is the momentum (in GeV/c) in the gap 
between hit k and k— 1; — 0.015a/ A^/X^., with the radiation length, 
is a constant depending on the composition and thickness of the tracking 
medium. Parametrizations for fj,^ are discussed further in Appendix O. The 
rms displacement in the length A& of 

S k = J\^, (4) 
V 3 p k 



is, for iron, given by 



10 GeV/c / A 



3/2 



<5 fe ~ 320 /im '- . (5) 

p V 10 cm/ 

For 10 cm tracking chamber separation, this displacement is the same as a 
typical spatial resolution measurement of a drift chamber. 

For constant chamber separation, one can set A k = A, = fx and 
incorporate energy loss effects in an approximate way to yield simplification: 

Sijip) ~ - A ^ — ^ | 2 min (i, j) 2 - 3 (i + j) min (i, j) + 6ij (6) 

^ / \ r 

H \~r) (min (i, j) + 1) 3min (i, j) — 4 (z + j) min (z, j) 

+i + j - min(i,j)]} , (7) 

with p the momentum at the start of the track. Many of the formulas 
presented here will assume the mean energy loss, (^/, is zero for simplic- 



ity, although, as will be seen, incorporation of finite (^^j can significantly 

improve momentum estimates from MCSQ 

1 For very low p tracks in long targets, muons can range out, in which case momentum 
determination from (gf ) is possible, in addition to the MCS estimate. 
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If the calorimeter is instrumented with a sufficiently large number of drift 
chambers, it is possible to exploit MCS to estimate muon track momentum 
from the scatter of the hits along a muon track. This can be seen from 
the following intuitive argument: Best estimates for 6 and yo follow from 
minimizing the x 2 '■ 

< yzi > - < yz >< ziz > 

"o — o > l°J 

< Z\Z\ > - < Z X ZQ > 2 

< z\Z\ >< yz > - < ziZq >< yz\ > 

V0 = ^7^^—^ • (9) 

Bracketed quantities < ab > are defined as 

aV- l b , N 

<ab>= — — - ; 10 
ZqV l z 

they are unchanged by a re-scaling of the error matrix. In the MCS limit, it 
follows that #o and yo are independent of the momentum, implying that x 2 °^ 
p 2 . Now, suppose one adjusts p until the x 2 probability density function 
attains its maximum. This occurs at x 2 — ^ 5 an d if the fit is performed at 
some nominal momentum po achieving x 2 = Xoi then an estimate for the 
true momentum of the track is 




P=P0\~2- (11) 



The MCS technique thus provides a momentum estimation method that 
does not require a magnetic field. 

2.2 Likelihood Function Method 

A more rigorous derivation begins with the observation that the joint proba- 
bility function for N correlated drift chamber hits can be written, assuming 
Gaussian errors, as 



1 \ N / 2 1 

P(y;0o,yo,p) = [tti n \ ex P 



7^x 2 {8o,yo,p) 



(12) 



2ttJ y/detV(p) 

with x 2 defined in equation |]. This can be converted to a (— )log likelihood 
function, 

£=±log(detV(p)) + ^ X 2 (Oo,yo,p), (13) 

where terms independent of y$, 6o, and p have been dropped. Estimates for 
yo, 9q, and p follow from minimizing C 
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- pV(po), where p 



2.2.1 MCS dominated limit 

In the multiple scattering limit, one can write V(p) 

is some nominal estimate of the momentum. From this, it follows that 

detV(p) = [pj detV(po) and X 2 (0o,yo,p) = ^X 2 (#o, Vo,Po)- In this 
limit, the log likelihood becomes, after dropping terms that are independent 
of 6q, y , and p 



C 



n 2 

P - 2/ 



-N\ogp + —^x (Oo,Vo,Po) 
z Po 



Minimizing with respect to the three fit parameters yields 



de 



8y 2 

dyo 



(14) 



0, as 



before, and equation |11. One can also obtain an estimate of the uncertainty 
in p from 

J_ = d^C _ 2 X 2 (8 ,y ,p ) 

r2 



a 



r d P 2 £=£ max Po 



(15) 



If the fit is iterated until p = Po an d the fit is reasonable so that \ 2 — N, 
then 



P 



1 



>2N 



(MCS limit). 



(16) 



2.2.2 Effects of Finite Spatial Resolution 

In the more typical case where chamber resolution is not negligible, one must 
solve the coupled equations 



d 2 



X (#0, 2/o, p) 



d_ 

Op 



d 

ot>o 

log(detV(p)) +x 2 {0 ,yo,p) 




0, 
0, 



(17) 
(18) 
(19) 



which can be accomplished via straightforward iterative methods by com- 
puter. 

Insight into the intrinsic resolution of the MCS momentum error estimate 
can be gained by examining an approximate expression for the momentum 
resolution, derived in Appendix [B]: 

-1/2 

\ (A 



p 



n=l 



2 \ 2 



\Sn + ) 



(20) 
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where £^ are the eigenvalues of the dimensionless scattering matrix S = 
^!i S(p).| For the geometry considered here, and ignoring energy loss, the 
relative momentum error is seen to be a universal function of the number of 
chambers n and the ratio p 2 /PmqSi 

^=F(n,p 2 /p 2 MCS ), (21) 

where 

PMCS = — , (22) 

defines a characteristic momentum scale (approximately 73 GeV for the 
NuTeV detector). Figure |l| shows plots of F(n,x) vs n for different values 
of x = p 2 /p\ics- About 7 chambers are required to measure pmcs to 50% 
fractional momentum error and 25 chambers to measure 10 x pmcs to 50%. 
For a given number of chambers n used on a track fit, one can define a 
critical value x cr it(n), such that a p /p < 50% for x < x cr i t (n). Figure || 
shows a plot of x cr it(n) vs n. For a given detector geometry, x cr i t (n) can 
be converted to p C rit( n ), the largest momentum that can be estimated from 
MCS scattering alone to 50% resolution. Figure ^ shows a plot of p C rit(n) 
vs n for the NuTeV detector. Momentum values of to 300 GeV can be 
estimated using 21 chambers in the detector, and up to 1 TeV using all 42 
chambers. 

The MCS limit is p 1 1 Pmcs ^ ^ or an n ' m wm ch case Eq. 16 is re- 



covered. If intrinsic chamber resolution dominates, MCS-based momentum 
estimates will then provide resolutions that behave as 

2 

(resolution limit),. (23) 



a v n~ ( P 



2N \pmcs 



2.3 Tracking in a Magnetic Field 

The much more straightforward way to determine momentum is via track 
displacement in a magnetic field. It is possible to combine the curvature 
measurement with the MCS momentum technique to improve the overall 
momentum estimation. 

For simplicity, the analysis will be restricted to a geometry of evenly 
spaced tracking chambers immersed in a uniform magnetic field oriented at 

2 This can be expressed in the computationally simpler form p 2 /(7p = 
2Tr[V(p)- 2 S 2 (p)]. 
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Number of Hits on Track 



Figure 1: Universal resolution function F(n,x) for MCS determination 
of momentum. The curves represent different values of x = which are 
x = 0.01 (solid- lower) , x = 0.1 (dots), x = 1 (dashed), x = 10 (dot-dashed), 
and x = 100 (solid-top). 
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8 




n 



Figure 3: Value of critical momentum P cr a, above which the fractional 
momentum resolution will exceed 50% for a given number of drift chambers 
n. This plot assumes the NuTeV detector geometry, with gq = 0.05 cm, 
A = 42.4 cm and 12.2 radiation lengths of material between each tracking 
chamber. 
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right angles to the track propagation. It is further assumed that the magnet 
Pt kick is much less that the momentum of the track being analyzed. In 
this case there is no dependence of the covariance matrix on fit parameters 
other than momentum, and the variance matrix for the fitted momentum 
takes the form 

/ 



E 



vv 



* 1 + 



1 



(24) 



where is given Eq. 
matrix, with 



2C and \l/ is the conventional spectrometer error 



-i 

12 
-1 
22 

-1 
13 

-1 
23 



33 



1 (p)z , 
i*iV _1 (p)z , 

k 
k 

A 2 - 

4p 



4 f 2 V 1 (p)z 2 , 



Here, 



(1,1,1,...,1) 

(^1,^2,^3, ••• 

( 2 2 2 
IZi, z 2 ,z 3 , .. 



zn) , 



, z N 



(25) 
(26) 
(27) 

(28) 
(29) 
(30) 

(31) 
(32) 
(33) 



and k = 0.0035, with B the magnetic field in Tesla assuming all spatial 
coordinates are in cm. 

In some detectors, such as NuTeV, the spectrometer follows the calorime- 
ter. Spectrometer momentum determination and MCS-based calorimeter 
determination are then independent and can be averaged. 

3 Results from Calculations 



3.1 Tracking in Unmagnetized Calorimeter 

Figures ||, [|, and |6| show results for the estimated fractional momentum 
error 5p = (a p /p) calculated for the NuTeV detector geometry from Eq. ^ 



as a function of various parameters appearing in Eq. 20 
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Figure |3| shows the dependence of fractional resolution on momentum for 
various numbers of drift chambers. Momentum dependence is present for all 
momenta and all numbers of chambers, indicating that the MCS resolution 
limit of <Tp/p = l/y2N is not reached until lower momentum. 



I 5 p (p k ,<J,D,X,14) 

Q . . • 

a- 

I 5 p [p k ,o,D,X,2lj 

I 5 p (p k ,o,D,X,28 

5 SJp, ,<J,D,X,35^ 

a _P\ k 

8 p (p k ,0,D,X,42^ 




1 10 100 1- 10 

Pfc 

Momentum (GeV/c) 

Figure 4: Momentum dependence of fractional momentum resolution for 
different numbers of chambers used in the fit. The curves correspond to 7 
chambers (upper-solid), 14 chambers (upper-dotted), 21 chambers (dashed), 
28 chambers (dot-dashed), 35 chambers (lower-solid), and the maximum 
possible 42 chambers (lower-dotted). This plot assumes the NuTeV detector 
geometry, with uq = 0.05 cm, A = 42.4 cm, and 12.2 radiation lengths of 
material between each tracking chamber. 

Figure |5| shows the dependence of cr p /p on the number of drift chamber 
hits. While the l/\/2N limit is not reached, the resolution does scale as 
A/y/2N, with A increasing with momentum. 

Figure ^ shows the dependence on chamber resolution. Effects are siz- 
able, indicating that a careful assessment of the intrinsic chamber resolution 
is necessary. 
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m 

Number of Chambers in Fit 



Figure 5: Dependence on number of drift chamber hits of fractional momen- 
tum resolution for different muon momenta. The dark solid curve represents 
the MCS limit. The other curves correspond to p = 5 GeV/c (lower-lighter 
solid), p = 10 GeV/c (lower-dotted), p = 20 GeV/c (dashed), p = 50 GeV/c 
(dot-dashed), p = 100 GeV/c (upper-solid), and p = 200 GeV/c (upper- 
dotted). This plot assumes the NuTeV detector geometry, with oq = 0.05 
cm, A = 42.4 cm, and 12.2 radiation lengths of material between each track- 
ing chamber. 
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0.55 



0.5 " 




Spatial Resolution (cm) 



Figure 6: Spatial resolution dependence of fractional momentum resolution 
for different muon momenta assumes 21 drift chamber hits in the NuTeV 
geometry. The curves correspond to p = 5 GeV/c (lower-solid), p = 10 
GeV/c (lower-dotted), p = 20 GeV/c (dashed), p = 50 GeV/c (dot-dashed), 
p = 100 GeV/c (upper-solid), and p = 200 GeV/c (upper-dotted). This plot 
otherwise assumes the NuTeV detector geometry, with A = 42.4 cm, and 
12.2 radiation lengths of material between each tracking chamber. 
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3.2 Tracking in Magnetized Calorimeter 

Figure |7[ H and ^ show the dependence of fractional momentum resolution 
in a magnetized calorimeter as a function of muon momentum, magnetic 
field, and number of chambers respectively. Also shown is the resolution 
estimate for a conventional momentum fit that incorporates MCS effects 
into the error matrix, but uses only the track curvature, not the pattern 
of scatter in the hits, to estimate momentum. The three plots assume a 
geometry with 0.05 cm resolution drift chambers separated by 10 radiation 
lengths of material. Figures |7| and ^ assume a magnetic field of 1 T, Figs. || 
and m assume a muon momentum of 50 GeV, and Figs. || and [|| assume 20 
chambers used in the fit. 



i 



»3 a pn fp, ,B,o,D,X,n' 1 

I r "\ k / 

I _ p k 

S a Pmcs (p k .B.o.D.X,N 




0.1 



10 



100 



M0 



Pk 

Momentum (GeV/c) 

Figure 7: Dependence of fractional momentum resolution on muon mo- 
mentum assuming twenty 0.05 cm resolution drift chamber hits spaced by 
10 radiation lengths of iron. The top (solid) curve assumes a conventional 
spectrometer fit, while the lower (dotted) curve incorporates information 
from MCS into the fit. 



Figure |10| shows the combined momentum resolution that can be achieved 
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Magnetic Field (T) 

Figure 8: Dependence of fractional momentum resolution on magnetic 
field assuming twenty 0.05 cmm resolution drift chamber hits spaced by 
10 radiation lengths of iron. The top (solid) curve assumes a conventional 
spectrometer fit, while the lower (dotted) curve incorporates information 
from MCS into the fit. 
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10 20 30 40 50 



n 

Number of Chambers 

Figure 9: Dependence of fractional momentum resolution on number of 
drift chamber hits for 50 GeV muon momentum and 0.05 cm resolution 
drift chambers spaced by 10 radiation lengths of iron. The top (solid) curve 
assumes a conventional spectrometer fit, while the lower (dotted) curve in- 
corporates information from MCS into the fit. 
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from the spectrometer and a varying number of calorimeter chambers used 
in the NuTeV experiment. The spectrometer alone provides a resolution of 
e s = 10%. 



o.i i 



i 



8„(p,.,a,D,X,7j e s 
1 

1 



1 



0.09 



1 + 1 

5„'p, ,o,T),X,2%) e„ 2 
P' k I * 0.085 



Muon Momentum (GeV/c) 



Figure 10: Dependence of fractional momentum resolution on muon mo- 
mentum for different number of chambers used in the MCS determination 
in a momentum estimate that combines the MCS estimate with a 10% spec- 
trometer measurement, es, in the NuTeV geometry (oo = 0.05 cm, A = 42.4 
cm, and 12.2 radiation lengths of material between each tracking chamber). 
The curves correspond to 28 chambers (dot-dashed), 14 chambers (dashed), 
and 7 chambers (dotted). The solid horizontal line at 0.1 represents the 
spectrometer-only resolution. 



4 Results of Monte Carlo Simulation 

The formulas developed in the previous section have been tested using a 
Geant simulation of the NuTeV detector (see Appendix |X|) using track find- 
ing and fitting algorithms described in Appendix This section presents 
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p,v, (GeV/c) 


20 


50 


100 


200 


P L view (GeV/c) 


20.6 


46.0 


93.4 


194.0 


P L ViewS (GeV/c) 


20.8 


46.5 


94.0 


195.0 


a/p l view 


0.078 


0.21 


0.30 


0.42 


cr/p^ v * ew (pred) 


0.17 


0.22 


0.26 


0.31 


; 2 views 
cj/p 


0.065 


0.15 


0.21 


0.30 


a/p 2 views (pred) 


0.12 


0.15 


0.18 


0.22 


e l p view (%) 


99 


99 


99 


96 


£ 2 views (%) 


99 


99 


99 


93 



Table 1: Momentum dependence of MCS momentum estimate from a Geant 
simulation of the NuTeV calorimeter. Results for reconstructed momentum 
(prec), resolution (a p /p), and reconstruction efficiency (e p ) are given assum- 
ing that either one or two views of drift chamber hits are used. The Monte 
Carlo sample contained 10 4 events, so statistical errors are of order 1%. 



results only for tracking in the unmagnetized NuTeV calorimeter.. 

Figures 11 and 12 show distributions of fitted values of 1/pasa function 
of track momentum using all 42 chambers in the NuTeV detector. Results 
are presented for momentum determination using only a single view in the 
drift chamber, and for fits that combine both views. Table summarizes 
momentum dependence of reconstructed momentum, fractional resolution, 
and tracking efficiency. 

Figures 13 and 14 show distributions of fitted values of 1/pasa function 
of the number of drift chambers used for a track momentum of 50 GeV/c. 
Results are presented for momentum determination using only a single view 
in the drift chamber, and for fits that combine both views. Table || sum- 
marizes chamber number dependence of fractional resolution and tracking 
efficiency. 

Agreement between observed resolution from the full reconstruction and 
Eq. ^ is satisfactory at 50 and 100 GeV/c. At 20 GeV/c, the observed 
resolution is considerably better than the prediction. Energy loss in the 
target is a significant fraction of the total muon energy in this case. The 
fitting procedure incorporates energy loss effects. Their inclusion introduces 
a second source of correlation between longitudinal chamber position and 
momentum that evidently enhances the resolution. Resolution at 200 GeV/c 
is about 30% worse in the Monte Carlo than predicted. The source for this 
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1/P-1 view, 42 chambers 



Figure 11: Distributions of 1/p (in (GeV/c) -1 ) estimated from MCS tech- 
nique for 20 GeV/c, 50 GeV/c, 100 GeV/c, and 200 GeV/c muons passing 
through 42 drift chambers in a Geant simulation of the NuTeV neutrino 
detector. Only one of two drift chamber views is used in the fitting. The 
histograms are for, from right to left, 20 GeV/c (solid), 50 GeV/c (dashed), 
100 GeV/c (dotted), and 200 GeV/c (dot-dashed) muons, respectively. The 
curves superimposed on the histograms represent simple Gaussian fits. 
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1 /P-2 views, 42 chambers 



Figure 12: Distributions of 1/p (in (GeV/c) -1 ) estimated from MCS tech- 
nique for 20 GeV/c, 50 GeV/c, 100 GeV/c, and 200 GeV/c muons passing 
through 42 drift chambers in a Geant simulation of the NuTeV neutrino 
detector. Both drift chamber views are used in the fitting. The histograms 
are for, from right to left, 20 GeV/c (solid), 50 GeV/c (dashed), 100 GeV/c 
(dotted), and 200 GeV/c (dot-dashed) muons, respectively. The curves su- 
perimposed on the histograms represent simple Gaussian fits. 
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1/P-1 view, 28 chambers 1/P— 1 view,42chambers 



Figure 13: Distributions of 1/p (in (GeV/c) -1 ) estimated from MCS tech- 
nique for 50 GeV/c muons passing through different numbers of drift cham- 
bers in a Geant simulation of the NuTeV neutrino detector. Only one of two 
drift chamber views is used in the fitting. The plots are for, clockwise from 
top-left, 7 chambers, 14 chambers, 42 chambers, and 28 chambers used in 
the fit. The superimposed curves represent simple Gaussian fits. 
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1/P-2 views, 7 chambers 1/P-2 views, 14 chambers 




0.06 



1/P-2 views, 28 chambers 1 /P-2 views, 42 chambers 



Figure 14: Distributions of 1/p (in (GeV/c) _1 ) estimated from MCS tech- 
nique for 50 GeV/c muons passing through different numbers of drift cham- 
bers in a Geant simulation of the NuTeV neutrino detector. Both drift 
chamber views are used in the fitting. The plots are for, clockwise from top- 
left, 7 chambers, 14 chambers, 42 chambers, and 28 chambers used in the 
fit. The curves superimposed on the histograms represent simple Gaussian 
fits. 
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N(chambers) 


7 


14 


28 


42 


/ 1 view 
cr/p 


0.56 


0.49 


0.30 


0.20 


/ 1 view/ i\ 
a/p (pred) 


0.52 


0.37 


0.27 


0.22 


/ 2 views 
cr/p 


0.44 


0.35 


0.22 


0.16 


a/p 2 views (pred) 


0.37 


0.26 


0.19 


0.15 


e 1 view (%) 


65 


91 


99 


97 


e 2 views (%) 


32 


82 


99 


95 



Table 2: Dependence of resolution, (— ) and reconstruction efficiency (e) 
on number of chambers used in the fit for 50 GeV/c input muons. Resulted 
are presented for fits using one or both views of the drift chamber. The 
Monte Carlo sample contained 10 4 events, so statistical errors are of order 
1%. 

disagreement is not fully understood, although it may be related to the 
significant tail that occurs on the high momentum side for fits to very high 
energy muons tracks. There are also small biases evident in the momentum 
reconstruction that, while much less than the momentum resolution, are not 
yet understood. For 50 GeV/c muons, the resolution is observed to scale with 
the number of chambers as 1/y/N, in agreement with the prediction. The 
expected y/T/2 improvement in resolution when combing the independent x 
and y views is also observed. 

5 Systematic Errors 

A proper survey of systematic errors requires treatment of real data, which 
will be presented in a forthcoming publication. A few obvious sources are 
commented upon here. 

The theory of multiple Coulomb scattering is well-established || [?], gj. 
The critical parameter that enters into fitting is /ifc, which represents the 
effective rms of a Gaussian approximation to the distribution of the projected 
scattering angle. As discussed in Appendix |C], different estimates for 
agree to within ~ 2% using simple parametrizations. It seems likely that 
this error could be reduced to negligible levels if a careful application of the 
Moliere theory is applied to a single material. 

Any drift chamber misalignment will produce a bias in the MCS fitting 
procedure that produces a fitted momentum estimate that is systematically 
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lower than the true value. Misplaced chambers will effectively introduce 
extra scatter between hits on a track. The only way the fitting routine can 
account for the extra scatter is to lower the momentum, thus increasing 
the contribution of MCS to V(p). If the misalignment is random, then the 
MCS fit will return relatively poor (high) values of the likelihood function. 
However, a correlated mis-alignment can mimic the effects of MCS fairly 
well. It is thus critical to have an accurate tracking chamber alignment. 

Spurious hits not directly associated with the muon track from chamber 
noise or multiple-pulsing of electronics will also produce undesirable scat- 
ter that will bias momentum fits towards too low values. Electromagnetic 
shower particles produced by high energy particles in dense calorimeters will 
produce similar effects. It is thus important to use only hits from "quiet" 
intervals of the muon track where there is no possible ambiguity. 

As is evident in Fig. || the resolution on the MCS fit is a fairly strong 
function of the chamber resolution a$. The MCS momentum estimate will 
also be biased by incorrect <to values in a positively correlated way, i.e., 
the fit will return too high values of momentum to reduce the MCS error 
contribution if ctq is input at a value larger than its true value. 



6 Conclusions 

It has been demonstrated that momentum estimates based on multiple 
Coulomb scattering can be extended to work for very high energy muons 
in dense calorimeters instrumented with a sufficient number of typical drift 
chamber tracking detectors. Using the MCS technique could allow high 
energy neutrino experiments such as NuTeV at Fermilab to increase their 
acceptance for low energy, wide angle muons that exit their detector before 
a spectrometer momentum measurement is possible. Other possible appli- 
cations exist in large detectors being assembled for long baseline neutrino 



oscillation searches, such as MINOS [14] at Fermilab. 

A The NuTeV Detector 

The NuTeV experiment will be used to provide specific examples of MCS 
momentum determination using analytic calculations and a detailed Monte 
Carlo simulation. NuTeV employs the Fermilab Lab E neutrino detector flTll, 



12] in a newly constructed high intensity sign-selected neutrino beam[ 



The experiment is currently taking data with the primary goal of measuring 



24 



the weak mixing angle sin 2 9w = 1 — Myy /M% to a factor of 2-3 times the 
precision achieved in previous neutrino experiments. NuTeV will also make 
improved measurements of nucleon structure functions and will perform high 
sensitivity searches for processes not predicted by the Standard Model of 
electroweak physics. 

The detector consists of a 690 ton iron target-calorimeter followed by an 
iron-core toroid spectrometer. The target-calorimeter is a 3 mx3 mxl8 m 
volume centered with its long side parallel to the neutrino beam. Hadron 
energy measurement, longitudinal event vertex determination, and event 
triggering and timing are performing using eighty- four 3 mx3 m x2.5 cm 
liquid scintillation counters read out through wavelength-shifter bars to pho- 
totubes on each of the four corners of the counter. The scintillation counters 
are separated by 10 cm of steel, and provide a sampling-dominated hadronic 
energy resolution of <je/E = 0.9/y/E + 0.3/E. 

Transverse event vertex position and muon track angle in charged current 
events are measured in the target-calorimeter using up to forty- two 3mx3m 
drift chambers which are spaced every 20 cm of steel. The drift chambers 
consist of two orthogonal planes, each containing twenty-four 12.7 cm wide 
cells. Each cell has two sense wires, permitting local left-right position 
ambiguity resolution. The chambers use a mixture of 50% argon- 50% ethane 
that produces a uniform 50 fim/ns drift velocity. Both sense wires of each cell 
are connected through shaping pre-amplifiers to multi-hit TDC's; the TDC's 
have 4 ns time buckets and can buffer up to 32 hits. When constructed 
and initially tested in the early eighties, the drift chamber resolution was 
measured to be 225 //m; the current resolution is estimated for studies in 
this paper to be do — 500 fim. Counting support structure, scintillation oil 
and other passive detector components, the drift chambers are separated by 
42.4 cm, corresponding to 12.2 radiation lengths. 

The toroid spectrometer the follows the calorimeter provides a 10% MCS 
momentum resolution for muons from 5 — 500 GeV/c. 

NuTeV detector response is simulated using a Geant-based package 0] 
that produces output in the same form as the on-line data acquisition system. 
Multiple scattering is simulated using the default Moliere scattering option 
of Geant. Muon energy loss in iron is simulated using the Landau fluctuation 
process for restricted energy loss. Important contributions from catastrophic 
energy loss from muon bremsstrahlung and e + e~ pair emission are carefully 
modelled by lowering Geant photon and electron energy tracking cuts to 0.1 
and 1.0 MeV, respectively. 
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B Derivation of MCS Momentum Resolution 



If one ignores energy loss effects, the track error matrix can be expressed as 

i 
p 2 



V(p) = a 2 I + ^S(^) 



If 



V(p ) \n) = \ n \n) , 
i.e., {A 2 } are the eigenvalues of V(po), then 



S(p )\n} = (\ 2 l -a 2 )\n) 



and thus, 



so that 



V(j>)|n) 



(a^-^+cx 2 



^0 /\2 
p 2 



log(detV(p)) = J2 lo & 

n 

Letting p = p + S p : 

log(detV(p)) ~ ^log(A 2 )- 2(A '' ^ 



AH 



PQ 



P0 



AS 



(34) 



(35) 
(36) 
(37) 
(38) 

(39) 



Expanding the x 2 function yields, for a linear fit, only one momentum 
dependence that remains after averaging: 



' ldx 2 (0o,yo,Po] 
2 dpdp 



5p5 p 




iozx - y z ) 



d 2 V(p) 
dpdp 



-i 



(40) 



x (y - Qqzi - VqZq)) SpS p 
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The trace can be expanded in the eigenvalues of V(po) : 



Tra 



dp 2 



-V(p) 



J2(n\^^— \rn) (rn\V( Pll )\n) (1.1) 



E 



2 (A« - ^) (A 2 - 4a 2 ) 



P 2 K 



Combining Eq. 39 with Eq. kill yields 



d 2 c _ (\i-giy 

d(5 P ) 2 V PoK 



(42) 



For a linear fit, there are no correlation terms between p and Oq or p and yoj 
so one has directly 



a 



(43) 



The inverse of the squared fractional momentum error is twice the sum of the 
ratios of squared eigenvalues of the MCS matrix to the total error matrix. 
One can alternatively write 



\2 2 _ A* c2 
P L 

where are the eigenvalues of the dimensionless reduced scattering ma- 
trix 



P 



S(p), 



yielding equation 20. 



C MCS Scattering Parameter 

The critical parameter in the error matrix is the quantity fj,k appearing in 
Eq. 44. The conventional PDG parametrization||, [K| 

K = 0.0136 J ^ (1 + 0.038 ln(A fc /X fc )) , (44) 

V A /c 

is the effective a of a Gaussian approximation to that part of the exact 
Moliere distribution that encompasses 98% of scattering angles. The NuTeV 
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drift chambers are separated by 12.2 ±0.1 radiation lengths, yielding a 
value of fif. of 0.0520 db 0.0002 GeV/c for single chamber separation. This 
value differs by (—0.7 + 3.41nn)% from the still-used simpler expression 
(jLk = 0.015-y/ Afc/Xo, where n is the separation between hits in numbers of 
chambers. The difference is small for hits separated by a single chamber, 
but becomes non-negligible for cases where multi-chamber gaps between hits 
develop due to inefficiency and noise. 



Lynch and Dahljl^] have analyzed the validity of Eq. 44 and give an 
alternate and more accurate parametrization in terms of the Moliere char- 
acteristic angle Xc and screening angle Xa, given by 



/;/ ^xc(k) = (4.0 x 10~ 4 ) ^Z k (Z k + l)p k A k /A k , (45) 



1 GeV/c 



rSvrc) Xaik) = ( 448 xl(r6 ) A 2/3 (! + 3-34^V), (46) 

where Z k and A k are the atomic number and atomic weight of the material 
in the scattering region, p k A k is the thickness of the region in g-cm 2 , p k 
is the momentum of the charge le ultra-relativistic particle, and a is the 
fine-structure constant. With these definitions, 



\ 2 X 2 M fl + v(k) 



p/(l GeV/c)J l + F 2 V v(k) 



In (1 + «(*:)) -l) , (47) 



where v(k) = 0.5O(A;)/(l — F), ft(k) = Xc(k)/Xa(k) is the mean number 
of scatters, and F is the fraction of scatters defined in a Gaussian fit to 
the Moliere distribution. Using Lynch and Dahl's prescription for including 
small effects of non-steel components of the NuTeV target yields, for F = 
0.98, fi k = 0.0528 GeV/c-for one chamber separation in NuTeV, 1.9% higher 
than the PDG formula prediction. 



D MCS Fitting Procedure 
D.l Hit Selection 

Initial track selection is performed using NuTeV's conventional calorimeter 
track finding and reconstruction software which has been tested over a period 
of many years in dozens of physics analyses. Inclusion of hits not associated 
with the muon track bias the MCS momentum fit by pulling the momen- 
tum to low values to handle the anomalously large scatter. These hits can 
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originate from the hadron shower induced by the neutrino interaction, from 
delta rays and other electromagnetic shower components produced by muon, 
and from electronic noise. To minimize their effect, only drift chamber cells 
with a single hit are used and a local trajectory requirement is imposed. 
The latter condition requires a hit position in a given drift chamber to be 
within 1 mm of the average position of hits from immediately neighboring 
chambers. One mm is about three times the typical single chamber MCS 
displacement for a 10 GeV/c muon. 



D.2 Energy Loss Correction 

A 50 GeV/c muon passing completely through the target calorimeter loses ~ 
17 GeV of energy on average. This degradation of energy causes downstream 
hits to spread further than would be predicted by the MCS error matrix 
calculated with a single momentum. This effect is taken into account by 
assuming an average energy loss of {dp/dz) = 0.0249 GeV/cm of detector, 
as calculated from the Geant simulation. No attempt is made to correct 
on an event-by-event basis for catastrophic energy loss, although such a 
procedure could be developed using the target scintillation counters. 

A smaller issue is the momentum p^ that should be used in Eq. ^. A 
simple analysis shows that this value is the geometric mean of the momentum 
evaluated at the beginning and end of scattering medium k 



Pk 



Po ~ \ -r 



(t) {Zk - l - Z0) 



(48) 



D.3 Fitting Algorithm 

Events with at least five drift chamber hits in at least one view are input 
to a likelihood minimization routine. The routine operates iteratively, vary- 
ing the momentum from iteration to iteration, but holding po fixed in the 
minimization of x 2 ($0) yo,Po) with respect to 9q an d yo- The momentum is 
set to 25 and 50 GeV/c for the first two iterations, and then minimized by 
approximating £ as a quadratic function of log p. A fit typically converges in 
8-10 iterations. A straight-through muon with 42 hits in each of two views 
requires about one second of CPU time on a 200 MHz Intel workstation. 

Figure [l5| show the likelihood function profile for typical fits to 20 GeV/c, 
50 GeV/c, 100 GeV/c, and 200 GeV/c momentum muons using 42 drift 
chambers. The shape of the function is approximately parabolic for a large 
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region about the minimum when plotted vs log p. Figure 16 shows C vs 
\ogp profiles for 50 GeV/c muons with 7, 14, 28, and 42 drift chambers used 
in the fit. As the number of chambers becomes small, the C distribution 
becomes skewed with a marked tail at high momentum. This behavior is 
expected since the relatively small and poorly determined scatter that occurs 
among a small number of chambers can be consistently described by nearly 
any momentum value higher than the true value. Note that C rises steeply 
for momentum values lower than the true value even for a small number 
of chambers. The MCS method may thus be used to set a useful lower 
bound on track momentum even if a reliable estimate of the true value is 
not feasible. 
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Figure 15: Typical likelihood function profiles plotted vs log p for (clockwise 
from upper left) 20 GeV/c, 50 GeV/c, 200 GeV/c, and 100 GeV/c muons 
fitted using 42 drift chambers in the NuTeV geometry. 
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Figure 16: Typical likelihood function profiles plotted vs logp for 50 GeV/c 
muons fitted using (clockwise from upper left) 7, 14, 42, and 28 drift cham- 
bers in the NuTeV geometry. 
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